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Abstract 

We investigate the problem of the existence of a density matrix P123 on a Hilbert 
space Hi &) H2 <8> % with given partial traces pu = Tr^ p\2z and P23 = Txi pnz- 
While we do not solve this problem completely we offer partial results in the form 
of some necessary and some sufficient conditions on pyi and p23- The quantum case 
differs markedly from the classical (commutative) case, where the obvious necessary 
compatibility condition suffices, namely, Tri pyi = Tt3P23- 

Mathematics subject classification numbers: 81V99, 82B10, 94A17 
Key Words: Density matrix, Entropy, Partial trace 

1 Introduction 

The problem considered here is closely related to the quantum marginal problem, on which 
there is an extensive literature. However, since the problem we consider involve overlapping 
marginals, results in the literature shed little light on it. We therefore introduce the 
problem in terms that we find natural, and postpone the discussion of the relation between 
our results and results on the quantum marginal problem until later in the introduction. 
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Let Hi, H2, H3 be three finite dimensional Hilbert spaces. Let pi 2 be a density matrix 
on H12 — Hi <g) H2 and, similarly, let p 23 be a density matrix on "H 23 = Hi ® r H-i- The 
question we ask, and which we can only partially resolve, is: 

Assuming that there is consistency of the partial traces, namely Tri p 12 = P2 = Tr 3 p 23; 
if/jai are necessary and sufficient conditions for the existence of a density matrix pi 23 on 
H123 = Hi <g> H 2 <8> H 3 swc/i i/iai Tri p i23 = p 23 and Tr 3 p i23 = p i2 ? 

This is an obvious question to ask in several contexts, e.g., [7]. We note the fact that 
in classical statistical mechanics the answer is that an extension always exists, and there 
is a simple formula for it: Given probability densities^ pi 2 (x,y) and p 23 {y,z) with 

Pi2(x,y)dfii(x) = J p2 3 (y,z)dp 3 (z) =: p 2 {y) , 

we may define 

/ \ Pi2(x,y)p 23 (y,z) , , 

Pi2 3 (x,y,z) = — . (1.1) 

p2{y) 

This extension is not unique, in general, but among all extensions, this one has the 
maximum entropy, namely Su + S23 — 5' 2 . The maximality is a consequence of the the 
Strong Subadditivity of Entropy (SSA) which says that for any extension 

5*12 + £23 > S123 + 5*2 • 

The classical SSA inequality is relatively easy to prove (as opposed to its quantum version); 
see [9]. The principle behind the construction is conditioning: Note that for random 
variables Y and Z with joint density p 23 {y,z), p 23 (F, z)/p 2 (y) is the conditional density 
for Z given Y . 

This construction of extensions by conditioning can be generalized to arbitrarily many 
factors. Given Pjj+i, j = 1, . . . , N, such that 



Pj-ij{x j - 1 ,x j )dfi j - 1 (x j - 1 )= I p j:j+1 (x j ,x j+1 )dfi j+1 (x j+1 ) =: pj(xj) 
define 

/ \ / \ TT Pj,j+l{ x ji x j+l) 
Pl,...,N[ X U ■ X N) = Pl,2{ X U x 2) 7 s • 

y PA x i) 

Again this extension is the maximum entropy extension by an interated application of SSA. 

However, the construction (11.11) . being based on conditional probabilities, does not gen- 
eralize to the quantum case, where the notion of "conditioning" has no obvious meaningful 
analog. In general, even for consistent density matrices, 

P12P2V23 



1 If the sample spaces over which x, y and z range are finite sets, these densities may be identified with 
diagonal density matrices in an natural way, embedding the discrete classical probability space problem 
into the quantum mechanical problem 
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need not be Hermitian, much less positive definite. While 

R = exp[logpi2 + logp 2 3 - logp 2 ] 

is Hermitian and positive definite its partial traces will not equal the desired reduced 
density matrices, even after renormalization, which would be required since the trace, 
Tfi23-R, is never greater than 1, and, generally, is less than 1. This fact follows from the 
triple Golden-Thompson inequality [8]: 



OO 1 1 poo 



Tri2 3 R < Tr 2 Tri 3 / p 12 — p 23 — d* = Tr 2 / p 2 — ^ — p 2 — ^ — d* = Tr 2 p 2 = 1 . 

JO t + p2 t + p 2 Jo t + p 2 t + p 2 

There are, in fact, consistent pairs of density matrices pi 2 and p 23 that have no exten- 
sion: Suppose that pi 2 is pure. If pi 2 3 is such that Tr 3 pi 23 = p 12 , the purity of pi 2 forces 
pi 23 to have the form p 12 ® p 3 and hence p 23 = p 2 ® p 3 . 

Thus, if p\ 2 is pure, the only compatible density matrices p 23 with which it has a 
common extension are the products p 2 <E> p 3 , in which case the unique common extension 
is P12 ® p 3 - Theorem 12. II in the next section generalizes this result by identifying a class of 
non-pure states p i2 , which can only be extended by product states pi 2 <8> p 3 . 

Motivated by these examples, and the obvious difference between the classical and 
quantum cases, we believe that a proper understanding of this problem will lead to a 
clearer understanding of entanglement and quantum information theory. 



1.1 The quantum marginal problem 

The quantum marginal problem, in its simplest form, is the following: Given density 
matrices pi and p 2 , on %i and % 2 respectively, let C(pi,p 2 ) denote the set of all density 
matrices pi 2 on Hi eg) 7-L 2 such that 

Trip i2 = p 2 and Tr 2 p i2 = p 1 . (1.2) 

C(pi,p 2 ) is the set of quantum couplings of pi and p 2 . 

Note that C(pi, p 2 ) is never empty; it always contains p\®p 2 . It is also evidently convex 
and compact, and hence it is the convex hull of its extreme points. Characterizations of 
the extreme points have been given by Parathasarathy [12] and Rudolph |14j . 

Other research has focused on the relations between the spectra of pi, p 2 and p 12 
that characterize the set of triples of density matrices satisfying (jl.2p . For results in this 
direction see the recent papers of Klyachko [6\, and Christandl et. al. [3], and references 
therein. 

The general quantum marginal problems concerns states on the product of arbitrarily 
many Hilbert spaces, and the marginals obtained by taking any combination of partial 
traces - for instance, the partial traces pi 2 and p 23 of pi 23 as in our problem. Most results 
pertain to the case of non-overlapping marginals, in contrast to the problem considered 
here. An exception is the paper by Osborne [TTJ. He uses the SSA of the von Neumann 
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entropy [10] to prove an upper bound the number of orthogonal pure states P123 such that 
Tri pi 2 3 = p 2 3 and Tr 3 pi 23 = pi 2 . We shall also employ entropy bounds, but are mainly 
concerned with the existence, or not, of mixed-state extensions P123 of compatible pairs P12 
and p 2 3- 

It is natural to use entropy in this investigation, and some early results on the two- 
space problem are given in terms of entropy. For example, if P12 is a pure state, then its 
partial traces pi and P2 have the same non-zero spectrum pQ. Conversely, if pi and P2 are 
two density matrices having the same non-zero spectrum, then there is a pure state P12 
satisfying (II. 2p . To construct it, assume without loss of generality that Hi and 7i 2 have 
the same dimension. There is a unitary U such that Up 2 U* = p\. Define \I> := Uy/~p~2, 
where the matrix on the right is regarded as a vector in "Hi ® H 2 . Then px2 = |^ / )(^ / | is 
such a pure state. 

If pi and P2 do not have the same spectrum, then the set C(pi,p 2 ) cannot contain any 
pure state. It is natural then to ask for the least entropy element of C(pi,p 2 ). Since the 
entropy is concave, the minimum will be attained at an extreme point, and the results of 
[I~2j [14"] are relevant, and lead easily to the minimizer in specific cases. However, a well- 
known inequality already provides a sharp a-priori lower bound for this minimum entropy 
coupling: 

The Araki-Lieb Triangle inequality p] says that when ( 11.2ft is satisfied, 

S 12 >\S 1 -S 2 \. (1.3) 

For an interesting discussion of entropy inequalities for a general classical marginal problem, 
see [I]. 



2 Necessary Conditions for the Existence of an Ex- 
tension 

Strong Subadditivity of Entropy (SSA) [TUl see also [9]] provides us with two necessary 
conditions for an extension to exist: Let P123 be an extension of p\ 2 and P23. Then, 
discarding the positive term S123 in SSA, 

Sl2 + S 2 3>S 2 . (2-1) 

Classically, one has monotoniciy of the entropy, meaning S\ 2 > S 2 and S23 > S 2 so that 
( 12.11) hold classically with 2S2 on the right side. However, quantum mechanically, equality 
can hold in ( 12. lft : Let P12 be a purification of P2, i.e., apure state |^ r )(^ r | on "Hi ®l-i 2 such 
that Tr 2 (|\l/)(\[ / |) = p\. Such purifications always existo Let P23 be the tensor product of 
p 2 and a pure state. Then S\ 2 = and S23 = S 2 . 

2 It is well known, and easy to see from the definition, that the set of all possible purifications ^ of pi 
is the set of all operators ^fp\V ', regarded as vectors in Hi <S> where U is a partial isometry from a 
subspace of H2 onto the range of pi . 



5 



A second form of SSA [10] leads to a sharper necessary condition: Again assume that 
P123 is an extension of P12 and P23. Then 

S12 + S 23 >Si + S 3 . (2.2) 

As we now explain, whenever ( 12.21) is satisfied by any consistent pair of density matrices 
P12 and P23 (not necessarily possessing a common extension), then ( 12.11) is automatically 
satisfied. To see this, note that by (II. 3p . S 2 — Si < Si 2 and S 2 — S 3 < S 23 . Adding these 
inequalities, we obtain 

25*2 < >Si2 + 5*23 + S\ + S3 . 

By (12. 2p . the right side is no greater than 2(Si 2 + S23), which yields (12. ip . 

The density matrices for which there is equality in the triangle inequality have a par- 
ticular structure: Let m and n be positive integers, and let {A 1; . . . , A m } and {pi, . . . , p n } 
be sets of positive numbers with Y^7=i^j = SlUi^fc = 1- Then, as shown in [2J, there 
exists a density matrix pi 2 such that the non-zero eigenvalues of pi 2 are {Ai, . . . , A m }, the 
non-zero eigenvalues of p 2 are {pi, . . . , p n } and the non-zero eigenvalues of pi are the num- 
bers {XjPk '■ l<j<rn,l<k<n}. For any such pi 2 , it is evident that S12 = Si — S 2 . 
Moreover, as shown in [2], whenever Si 2 = Si — S 2 , the spectra of pi 2 , pi and p 2 are related 
in this way. 

The following theorem generalizes the observation that pure states pi 2 may only be 
extended by product states pi 2 (g> p 3 . When P12 is pure, — Si 2 — 5*2 — Si. 

2.1 THEOREM. Let p 12 be a density matrix such that 

Si 2 = Si — S 2 . (2-3) 

Then p 12 and p 23 have a common extension if and only if p 23 = p 2 (g) p 3 . 

Proof. Using (T23D in ([OD we obtain 

Si — S 2 + 5*23 > 5*3 + Si . 

which reduces to S23 > S 2 + S3. By the subadditivity of the entropy, this means that 
5*23 = S 2 + 5*3, and so p 23 = p 2 <g> p 3 . □ 

In Section 4 we give an example in which (I2.2p is satisfied, but there is no common 
extension. 

3 Sufficient Conditions for the Existence of an Exten- 
sion 

We do not have any very general sufficient conditions for compatible pairs to possess a 
common extension. One general positive statement that can be made is the following: 
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3.1 THEOREM. Let pi 2 and p 23 be a compatible pair of density matrices that posses a 
common extension pi 23 that is positive definite. Let || • || denote the trace norm. Then there 
is an e > 0, depending on the dimensions and the smallest eigenvalue o/pi 23; such that if 
P12 and p 23 is another compatible pair on the same spaces and 

11/012 - P12II + H/O23 - P23 1 1 < e , (3.1) 
then P12 and p 23 possess common extension. 
Proof. Define 

Pl23 = Pl23 + [p~12 - Pn] ® P~3 + Pi ® [p 2 3 - P23] + Pi ® [p2 ~ P2) ® P~3 • 

It is easily checked that Tr!p 12 3 = p 23 and Tr 3 p 123 = p~i 2 . Furthermore, p~i 23 is self-adjoint, 
and under the condition ( 13. II) . is positive when e is chosen sufficiently small. □ 

3.2 Remark. Of course, in this finite dimensional setting, the trace norm could be replaced 
by any other norm. Also, Suppose that pi 2 is positive definite, and p 23 = p 2 <E> p 3 with 
p 3 positive definite. In this case p i2 and p 23 are compatible and have the positive-definite 
common extension p 12 <g> p 3 . The theorem says that any compatible pair that is a small 
perturbation of pi 2 and p 23 has a common extension. However, as Theorem 12.11 shows, the 
requirement of positive definiteness cannot be dropped. 

Let p u and p 23 be a consistent pair of density matrices, both of which are separable. 
Recall that a bipartite density matrix p 12 is finitely separable if is is a convex combination 
of product states; i.e., if 

N 

P12 = £ \A 3) ® # 

where each cr^ is a density matrix on each t% is a density matrix on "H 2 , and, 
crucially, each Xj > 0. The set of separable density matrices is the closure of the set of 
finitely separable density matrices. 

Separable density matrices are often viewed as being a classical ensemble of product 
states. As such, separable states are free of quantum correlations among observables on Hi 
and H,2- Separable bipartite density matrices do, indeed, behave much more like classical 
joint probability distributions. For example, if p i2 is separable, then (see e.g. [2]) 

S\2 > max{S'i S 2 } , 

and thus when both p i2 and p 23 are separable and compatible, S12 + 5* 23 > Si + S3, and 
thus the condition (I2.2p . which is necessary for a common extension to exists, is always 
satisfied. 

One might therefore, hope that a common extension p 123 would exist whenever p 12 and 
p 23 are compatible and both are separable. In the next section, we show that this is not the 
case: Further conditions must be imposed to ensure the existence of an extension. Here is 
one case in which a common extension does exist: 
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Suppose that 

n n 

p 12 = Aj P (i) ® and p 23 = ^ (8) r y) (3.2) 

i=i i=i 

where the and r^-* are density matrices and the Aj are positive numbers with 
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Pi23 := 2^ A i ^ ® (T °' ) ® r(i) ( 3 - 3 ) 

j=l 

is a common extension. 

While this example is based on a strong assumption, namely that the weights and the 
factors on % 2 coincide, we shall show in Section 4 that consistency and separability is not 
enough. 

However, as we now explain, there is a more general version of this construction using 
coherent states. Good references on coherent states are [5l [13] . 

Let H be any d dimensional Hilbert space. Define J = (d — l)/2. Then there is an 
irreducible representation of 577(2) on TL. Associated to this representation is a family of 
pure states |0) (0| on %, parameterized by a point Q in the unit sphere, S 2 . Given any 
self-adjoint operator A on H, there is a function a(Q) on S 2 such that 



A= / dfi a(fi)|n)<n| (3.4) 
Js 2 

where the integration is with respect to the uniform probability measure on S 2 . The 
function a(O) is unique if we further require that a(Q) be a spherical harmonic of degree 
no higher than 2 J. (See [5J Page 33 and Equation (4.10)].) This function a(Q) is called 
the upper symbol of the operator A. The upper symbol of A need not be a non-negative 
function on S 2 even if A is positive definite. However, if a is non- negative, then the operator 
A defined by ( 13 .4p is positive semidefinite. 

Now consider a bipartite density matrix p 12 on Hi ®H 2 where the dimension of Hj = dj 
for j = 1,2. Suppose that has a representation 



9x2= \ dn 2 dfi 1 p 12 (fi 1) n 2 ) (ifliXfiii® \n 2 )(n 2 \) 

Js 2 Js 2 

with ^12(^1, ^2) a nonnegative function on S 2 x S 2 that is a spherical harmonic of degree 
dj — 1 in Qj, j = 1, 2. 

Then taking the trace over Hi, we find 



p 2 := Trxp 12 = / dQ 2 
Js 2 



dQi pi 2 (tti,Q 2 ) 
Us 2 



\n 2 )(n 



2| ) 



and evidently / dOi pi 2 (£li, ^2) has the degree d 2 — 1 in tt 2 , which ensures that it is the 
Js 2 



's 2 

unique upper symbol of P2 of minimal degree. 
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Likewise, if P23 is a density matrix on r H. 2 <8> "%3 with the dimension of 7i 3 = dj, and if 
P23 is compatible with P12 and moreover 



p 23 = / d^ 2 / dn 3 p 23 (n 2 ,n 3 ) (\n 2 )W®\£i 3 )(n 3 \) 

Js 2 Js 2 

with ^23(^2,^3) being a nonnegative function on S 2 x S* 2 that is a spherical harmonic 
of degree dj — 1 in Qj, j = 2,3, f s2 dfl 3 p 23 (Q 2 , ^3) is the upper symbol of p 2 , and thus 
quantum compatibility of P12 and P23 implies the classical compatibility of the probability 
densities pi 2 (^i, fl 2 ) and p~23(^2, f2 3 ). Hence the classical prescription may be used to define 

~ (n n n \ £12(^1, ^2^23(^2,^3) 

Pl23l"l,"2,^3) := 



p 2 (^ 

and it is then the case that 



P123 



dfti / an 2 [ dQ 3 Pi23(Oi,Q 2 ,o 3 ) (iniXHii ® |n 2 ><02[ ® |n 3 ><fi 3 |) 

Js 2 Js 2 



>s 2 Js 2 Js 2 
is positive semi definite, and is a common extension of P12 and P23. 



4 Separable Compatible Pairs with no Extension 

In this section we give examples of compatible pairs P12 and P23 that do not have an exten- 
sion, but nevertheless satisfy the necessary condition (I2.2p . Moreover, in these examples 
both P12 and P23 will be separable. As we have remarked above, one might expect the 
extension problem to simplify in the presence of separability This is not the the 
following examples show. For simplicity, our examples will be on C 2 <E> C 2 <8> C 2 , but of 
course may be embedded into higher dimensional spaces. 

4.1 LEMMA. Let p 2 be a non-pure density matrix on C 2 ; i.e., the rank of p 2 is 2. Let 

{•01, ^2} be an orthonormal basis of C 2 consisting of eigenvectors of p 2 . Let p 2 ipj = Pjipj, 
j = 1,2. Then there exist two unit vectors <j>\ and <ft 2 in C 2 and positive numbers v\ and 
v 2 such that 

2 2 

i=i 3=1 
and such that the four vectors ip 2 , and <p 2 are pairwise linearly independent. 

Proof. If P2 is |l, we may take {tpi, ip 2 } to be any orthonormal basis of C 2 , and then choose 
02 } to be any other orthonormal basis such that (f)\ is not proportional to either ip\ or 
ip 2 . Then the four vectors <pi, cj) 2 , ifti, ip 2 are pairwise linearly independent. In this case, we 
take pj — Uj — I for j = 1,2. 

Next, assume P2 has distinct eigenvalues p\ > p 2 . Let {ipi, tp 2 } be an orthonormal basis 
consisting of eigenvectors of p 2 . Let <f>\ be any unit vector that is not proportional to either 
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ipi or i[) 2 . Then there is a unique largest number v\ > so that p 2 — ^i|</>i)(0i| is positive 
semidefinite. Thus, p 2 — vi\(j>i) (4>i\ is rank one, and hence 

P2 ~ ^l|0l)(0l| = ^2|02)(02| 

for some unit vector <ft 2 that is not proportional to <j>\ since p 2 has rank 2. Likewise, 2 
is not an eigenvector of p 2 , since by the non-degeneracy, this would force to be an 
eigenvector, which it is not. □ 

4.2 THEOREM. Let p 2 be a non-pure density matrix on C 2 . Then there exist separable 
density matrices pxi and p 23 , extending p 2 , and such that (12. 2p is satisfied (with strict 
inequality) , but such that p\ 2 and p 23 have no common extension px 23 . 

Proof. Let {tpi, ip 2 } and {</>i,</> 2 } be sets of unit vectors such that (11. 2ft is satisfied, and 
where {ipi,ip 2 } is an orthonormal basis of C 2 consisting of eigenvectors of p. 

Let {xi, X2} be any orthonormal basis of C 2 , and let {771, rj 2 } be any linearly independent 
set of unit vectors in C 2 . Define 

2 2 

P12 = ^2^j\Vj ® 1>j)(Vj ® V'il and P23 = "jI^j ® ® Xj\ ■ (4-2) 

i=i i=i 

Then evidently Tri p i2 = Tr 3 p 23 = p 2 . Next, for any vector (z,w) G C 2 , let (z, w) ± = 
(—w,z) so that (z,w) is orthogonal to (w,z). 

Suppose a common extension p 123 does exist. Then rj^ <g> ^ is in the nullspace of pi 2 , 
and hence 

2 

5^(f/T ® V>i ® Xj, P123 Vi®^® Xj) = (Vi ® ^1, P12 »7i" ® ^1) = . 
Therefore, since pi 23 is positive, both 77^ <g) ^1 <S) Xi and 77^ <g> ipi g) X2 are in the nullspace 

Of P123- 

Since r) 2 ®^ 2 is in the null space of pi 2 , the same argument shows that both r] 2 ®i}) 2 ®Xi 
and rj 2 ®ip 2 ® X2 are in the nullspace of pi 23 . Thus, the four vectors 

r]i Xj and r] 2 <g> ^2 <8> Xi > j = 1, 2 

belong to the nullspace of p i23 . 

Likewise, both (pf- ® xi and 4> 2 ® X2 belong to the nullspace of p 23 . Arguing as above, 
we see that the four vectors 

Tjj- ® ^i 1 ® xi and ^ ± ®0 2 L ®X2, J = 1, 2 

belong to the nullspace of pi 23 . 
Define the unit vectors 

#1 = ^^i^Xi , ^2 = Vi®i J i®X2 , ^3 = vi~ < 8>i)2<8>Xi > and * 4 = r) 2 <S>ip 2 <S>X2 , 

(4.3) 
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and the vectors 

$1 = Vi®4>i®Xi , $2 = 7^802 ®X2 , $ 3 = V2®4>i®Xi , and $ 4 = V2®4>2^X2 , 

(4.4) 

These 8 vectors are in the nullspace of plunder the hypotheses imposed so far. 

Now let us temporarily impose the additional hypothesis that {771,772} is orthonormal. 
Consequently, {^i - ,^} is orthonormal. 

We claim that the 8 unit vectors *ffj, $j, j = 1, 2, 3, 4, are linearly independent. To see 
this, note that under our additional hypothesis, {^1, ^2, ^3, ^4} and {$1, $2, $3, $4} are 
orthonormal, and $j) = 8^. Therefore, if 

4 4 

j'=i i=i 
(oj^j + 6j-$j) = for each j = 1, 2, 3, 4. 

However, since the vectors ipi and are linearly independent, so are the vectors ip± and 
<t>2 = 0i ■ Thus ^1 and $1 are linearly independent, and hence 01 = 61 = 0. 

The same sort of argument shows that Oj = bj = for each j = 1,2,3,4. Thus the 
eight vectors \ET,-, j = 1,2,3,4 and $j, j = 1,2,3,4 are linearly independent and in the 
nullspace of p i2 3, Since the dimension of C 2 <S> C 2 <S> C 2 is 8, this means pi 23 = 0, which is 
impossible. 

Therefore, the compatible pair of density matrices pi 2 and P23 defined by (14. 2 j) has no 
common extension. 

Next, observe from (I4.2[) that 

2 

i=i 

Since 771 and r]2 are orthogonal, and since r]i <8> ^1 an d ^2 <E> ^2 are orthogonal 

2 

5*1 = S12 = -J^pj logpj . 

i=i 

In the same manner we see that 

2 

5*3 = S23 = - ^2 v i log v i ■ 

Thus we have 

5*12 + 5*23 = Si + S3 . 

Finally, making a sufficiently small change in {771,772}, we may arrange that this set 
of unit vectors is no longer orthogonal, but that the 8 unit vectors *f?j,$j, j = 1,2,3,4 
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defined in (14.31) and (14. 4p are still linearly independent. Then the compatible pair pi2 and 
P23 defined by (14.21) has no common extension. It is still true that S23 = S3. Also, since 
{^1, ^2} is orthonormal, {r]\ <g> ipi, r) 2 Cg) ip 2 } is orthonormal, and so S\ 2 = — Y^=i l°g 
However, pi = Y^j=\ ^j\Vj){Vj\ and since {771,772} is not orthogonal, 

2 

St < -^fijlogfij = S 12 ■ 

Therefore S12 + 523 > + S 2 . □ 

4.3 Remark. While the condition that p\ and p 2 have the same non-zero spectrum is 
necessary and sufficient to ensure that there exists a common purification of p\ and p 2 , the 
above construction shows that the condition that the nonzero spectrum of pi 2 equals that 
of P3, and that the nonzero spectrum of P23 equals that of pi, together with compatibility, 
does not ensure that there exists a common purification of p\ 2 and p 23 . Indeed, let {r/i, r/ 2 }, 
{^1, 1P2}, {0i, 02} and {xi, X2} be four orthonormal bases of C 2 having no vectors (or their 
opposites) in common. Then with pj — i/j — 1/2, j — 1, 2, (14. 2p defines a compatible pair of 
density matrices with p 2 = The non-zero spectrum of pi 2 , P23, Pi and p 3 is {1/2, 1/2}, 
and yet there is no common extension of P12 and P23, as the proof of Theorem 14.21 shows. 
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